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We study the production of a pseudoscalar Higgs boson A^ in association with a Z" boson at a 
future international linear colhder (ILC). We consider the contributions to this process at the one 
loop level in the Minimal Supersymmetric Standard Model (MSSM) from top and bottom quarks 
as well as stop and sbottom squarks. We also study the squark contributions to the decay widths of 
the pseudoscalar Higgs boson for the decays A" — > 'yZ° and A*^ — > Z°Z°. The contribution from the 
supersymmetric loops are found to be directly proportional to the squark mixing and potentially 
If^ ' large due to the massive pseudoscalar Higgs coupling to squarks. 
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^ ■ I. INTRODUCTION 

<N 

' The Higgs mechanism is the means by which the electroweak symmetry is broken in the Standard Model (SM) 
^ '. and in the Minimally Supersymmetric Standard Model (MSSM) (see [111 11 for review). The MSSM has two Higgs 
^ • doublets which are used to generate the masses of the up- and down-type quarks. This leads to five physical Higgs 
[ bosons. These consist of two CP-cvcn neutral scalar bosons (/i°, 77°). one CP-odd neutral pseudoscalar (^°), and 
• two charged bosons (H^). The extended Higgs sector has several new parameters that can be determined with the 
' input of two parameters, the mass of the pseudoscalar M^o and tan/3 given reasonable assumptions as to the size 
of the other supersymmetric parameters in the theory. The parameter tan/3 = Vu/vd is the ratio of the vacuum 
expectation values (vevs) of the up and down sectors. We are interested in the phenomenology of the pseudoscalar 
Higgs boson at a future international linear collider, in particular, what the heavy squark contributions can tell us 
about the phenomenology. 

[ There have been many studies of pseudoscalar Higgs boson phcnomcnology0,0,0,ll,0E,0|. The phenomenology 
of the pseudoscalar Higgs produced in association with a Z° boson at a hadron collider (such as the Tcvatron 
Q^; and the CERN LHC) has also been studiedUl 111 EHIl El, as well as at a future international linear collider 
D . (ILC) \rA IT9I I20I I2H . There have been some studies of the decays of Higgs bosons into squarks p^. |2^ that show 
' how large these processes can become, but squark contributions are rarely taken into account for pseudoscalar Higgs 
processes. In this paper we will discuss squark contributions to processes involving the production and decay of a 
pseudoscalar Higgs boson which have been missing from the literature thus far. 
. , The construction of an ILC would greatly enhance our ability to measure the parameters in many new physics 
' scenarios 13 m including the MSSM. The process e+e~ Z^A^ at an ILC is interesting because it has the possibility 
of four bottom quarks in the final state allowing for precise measurements of the process as well as many other 
excellent final states. This process is also interesting because it could compete with the e~^e~ — > process. The 

e+e~ — > Z^A^ process occurs at the one loop level because the pseudoscalar Higgs boson does not couple to vector 
bosons at tree level. This process also allows for the exploration of squark mixing for the heaviest two generations of 
squarks. 

We find that the dominant contributions to the e+e^ Z^ process does not come from top and bottom quark 
loops. We find that the squark contributions to the e+e^ Z^ process are relevant at all values of tan/3 and 
although they depend on the mixing in the stop and sbottom squark sectors, they dominate over the standard model 
field contributions due to the large coupling of the pseudoscalar to squarks. We find that the pseudoscalar decay 
— > Z'^Z'^ has significant contributions from the squark sector and that the — > Z^Z^ branching ratio becomes 
more important for large values of the mass of the pseudoscalar. We also find that the jZ'^ decay squark has 

contributions on the order of the quark contributions for a light pseudoscalar Higgs at large values of tan/3. Overall, 
squark contributions need to be added to processes involving neutral Higgs bosons to complete our understanding of 
their phenomenology. 
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II. SQUARK CONTRIBUTIONS 

111 the SM, quarks generically couple to the pseudoscalar Higgs boson as ^^mq/v. Beyond this base couphng, the 
up- and down-type quarks couple to the pseudoscalar differently. Up-type quarks couple as cot/3 and down-type 
quarks couple as tan/3, leading to the well-known conclusion that bottom quarks become more important at large 
values of tan/3. This is also true in the stop/sbottom sectors. The sbottom squarks become more important as tan/3 
becomes large. 

In the MSSM, right and left handed quarks, Qr^l, have scalar super-partners, qh.R- We are interested in the stop 
and sbottom squarks (^ij, ii, fo_R, The squark sector of the MSSM has the possibility for squarks to mix into mass 
eigenstates that arc different than the usual left-right basis. We can introduce a mixing angle that diagonalizes the 
squarks into their mass eigenstates. This can be written generically for the new squark eigenstates gi^2 as0 

qi= QLCq+qFlSq, 

q2 ^ -QLSq + QRCq, (1) 

where Cq = cos 6q and Sq = sinOq are the mixing angles in each of the squark sectors. Squark mixing is particularly 
interesting when it comes to the squark loop contributions for pseudoscalar Higgs production. In the left-right basis, 
the pseudoscalar Higgs couples only to squarks that change their handedness at the vertex leaving only vertices such 
as A^qRqi^. However it is now easy to see that in the presence of squark mixing, vertices of the type A^q^qj for all 
values of («, j) are present and retain the size of the A^qRq^ vertex times a function of the new mixing angle. This is 
also true of other neutral Higgs bosons in the MSSM, but the analysis is more complicated as there are more vertices 
present for the and Higgs bosons. 

In this model, all of the pseudoscalar Higgs vertices in the squark sector are proportional to the squark mixing 
angles. Thus the mixing angles regulate the squark contribution to any pseudoscalar Higgs process. It is possible to 
tune the mixing angles Sq and Cq to eliminate all but one of these vertices. For instance if sin^t = 0, then only the 
A^tit2 vertex remains and it is maximal. It is easy to see in Eqn. that this would correspond to being back in the 
left-right basis. 

The mixing angles themselves are not completely independent from the other parameters in the MSSM. We can 
write the mixing angles as 

^ 2mtop(^-Mcot/3) 
mi - mi 

El E2 

• (^Q\ 2TObot(A - Mtan/3) 
sm(26';,) = 2 2 ■ (3) 

TO- - TO- 

Here At^b are the tri-linear scalar couplings for soft supersymmetry breaking in the MSSM Lagrangian and fi is the 
supcrsymmetry breaking Higgs mass parameter in the super potential. 

We can see that there is more splitting in the stop sector as it is proportional to the top quark mass which is much 
larger than the bottom quark mass. This is a general feature of squark mixing. In our analysis we picked toi = 1 TeV 
in both the stop and sbottom sectors and generated TO2 using the other parameters in the theory for an arbitrary 
mixing angle allowing us to study the mixing in the squark sector. Thus we always have TO2 > mi in the stop sector. 
This is the standard ordering for squark masses. Our choices for squark masses are well beyond current experimental 
limits |2^ and we have chosen such heavy squarks to make our results more conservative. The other parameters in 
the theory were chosen to be ^ = 300 GeV and At^b = 1500 GeV. Of course, the squark contributions to the process 
depend on the mass of the squarks and lighter squarks lead to larger contributions to any process involving them. 
In choosing fhi = 1 TeV we find maximal mixing in the squark sector leads to the lightest masses for the second 
quark. Although, our results show a very complicated dependence on these mixing angles, the dominant effect of the 
supersymmetric contributions are seen when the squarks are the lightest as is to be expected. 

These mass relations are tree level relations and the corrections to the mixing angles are known in the literature |27j. 
However, since we are interested in studying the effects of the squark mixing on our observablcs, it is unclear how to 
incorporate these corrections to produce squark masses that obey specific mixing angles. The effects of these higher 
order corrections are the most pronounced in the sbottom sector and future studies should be completed on this 
subject. We have used the MS running-mass of the bottom quark in our generation of the second sbottom squark 
mass. 
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TABLE I: Feynman rules for squark mass eigenstate couplings to the photon, a Z boson, two Z bosons, and the pseudoscalar 
Higgs boson. Here the squarks have momentum p and q (running in the direction of the flow of charge) and the photon and 
Z" boson have a free Lorentz index fi. The two Z" bosons have free Lorentz indices {fi, v). The (i, j) indices are written in the 
direction of the flow of charge at the vertex. is the isospin of the left-handed quark for which the squark is the super-partner, 
±i, Qq is the charge of the squark, = sin6'uj, and Cw = cosO^. The photon couples diagonally to the mixed squarks. The 
pseudoscalar coupling constant is listed in the text and is potentially very large compared to the coupling in the quark sector. 

III. PSEUDOSCALAR COUPLING TO SQUARKS 

Squark mixing imposes new Feynman rules for the vertices in our processes. The results of the apphcation of the 
squark mixing to the pseudoscalar vertices can be found in Tabled The coupling of the squarks to the pseudoscalar 
can be written as 

^ = -^(/i-A4Cot/3), (4) 

V 

i;, = -^(M-Atan/3), (5) 

V 

where v = + w| = 246 GeV is the Higgs VEV and we have chosen sign(/i) < implicitly in our choice of couplings. 

The squark contributions to pseudoscalar Higgs processes arc often missing in the literature and were believed to 
be quite small because they are loop suppressed and the squark would be very heavy if they exist in nature. Although 
this is certainly a possibility given the right parameter choice, the pseudoscalar coupling to squarks is quite large and 
is not entirely suppressed by the heavy squark loops. This is due to the fact that the pseudoscalar couples to right- 
and left-handed squarks only in a non-diagonal way. In the MSSM, the non-diagonal squark vertices with any neutral 
Higgs boson ($ = {hP , H'^ , ^"^Dj written here as ^gnQL, is proportional to the tri-linear soft-supersymmetry breaking 
parameter multiplied by the same factor as the coupling in the quark sector modulo any mixing angles that may also 
be present at the vertex. For the pseudoscalar Higgs, there are no mixing angles present, just tan/3 or cot/3 depending 
on the nature of the squark. This could also increase this effect in the sbottom sector once the difference in the 
masses between the top and bottom quarks is overcome. Therefore the squark-squark-Higgs coupling is potentially 
greatly enhanced compared to the quark-quark-Higgs coupling, so much so as to overcome the loop suppression and 
the heavy squark masses to become the dominant contribution. 

If we consider the top/stop sector as an illustration, the importance of the squark contributions will become clear. 
The vertex A'^tt ~ mt / v and the dominant contribution to the vertex A^tt ^ rrit /v x At cot /3. Thus, the squark vertex 
is as large at the quark vertex times a large factor, in our case 1500 GeV for tan/3 = 1. This factor of 10^ is squared 
in a cross-section or width and gives an overall factor of 10^. We have picked squarks with masses on the order of a 
TeV (about five times the mass of the top quark) which deceases the cross-section or width by an order of magnitude 
when compared to a squark with a mass lowered to that of the top quark. So naively we would consider the squark 
contributions to be larger than the quark contributions by a factor of 10^. This is exactly what wc have found for the 
e+e^ — > Z^A'^ process and the A" — > Z^Z" width. There are some cancellations in the A" width leading to a 

result that is on the same order as the contributions from the SM fields for some values of Af^o and tan /3. 

We would like to emphasize that this is the case for the non-diagonal squark vertices in the entire neutral Higgs 
sector of the MSSM. The hPqB,qL and H^q^qL should show similar enhancements with some additional dependence 
on the mixing angles (a,/3). The squark loop processes involving the production of these particles could very well 
be larger than the tree-level production processes and could lead to better bounds on squark masses or trilinear 
soft-supersymmetry breaking terms from the LEP2 data for the production of the other neutral Higgs bosons in the 
MSSM. 

It should be noted that this is partially an arbitrary enhancement because we have chosen the soft-supersymmetry 
breaking parameters to be large. However, if the soft breaking terms exist in nature, they will enhance the squark- 
squark-Higgs vertex over that of the quark-quark-Higgs vertex proportionately to their size. With this information we 
are ready to construct the matrix elements for our decays A" 7Z" and A'^ Z'^Z^ and our process e+e" Z^A^. 
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FIG. 1: The three-point functions needed for the e^e^ —^ Z"A calculation. Top and bottom quarks as well as stop and 
sbottoin squarks were allowed to run inside the blob. The momentum assignment was chosen to make the e^e~ Z'^A'^ 
calculation simpler. 



IV. MATRIX ELEMENTS AND RESULTS 



To understand the e+e~ — > Z'^A'^ process, we calculated the three point functions jZ'^A^ and Z'^Z^ Al' with both 
top and bottom quarks and stop and sbottom squarks. A generic representation of the three-point functions can be 
seen in Fig. (^). The contributions from the top and bottom quarks were checked against known results|23. We also 
checked our calculation with the known partial width for Vi{Ai^ — > 7^°) as given by hdecay|23. We found excellent 
agreement. We also found that due to the tensor structure of the quark and squark contributions, the two processes 
do not interfere with each other. The exact form of the contributions is worked out in the appendix and it is shown 
explicitly that the two contributions do not mix. 

The squark contributions to ^ii^J^ 7^") become less and less important as M^o becomes large. This can be 
easily seen in our parameterization of the matrix elements squared below 



\M 



SUSY I 2 

A° 



1 



167r2 (167r2)2 



(M2 



AO 



Ml? 



(6) 



where the superscript SUSY here implies that this is the contribution from the squark fields (not the SM fields) and 
Nc = 3 is the number of colors. 

The functional form of the i?* functions can be found in the appendix. However, it is simple to see in Eqn. © 
that the second term makes the amplitude smaller and smaller as M^o ^ Mz- The functional form of the A'i^^E'i^ 
functions guarantee the matrix elements are positive definite. What cannot be immediately seen is that the two 
pieces of the squark contribution almost cancel. Even though the pseudoscalar squark vertices are greatly enhanced 
compared to their standard model counterparts due to the trilinear term, the supersymmetric contribution to the 
width is only slightly larger than the SM contribution in some of the parameter space. 

We also needed to calculate the three point function for Z^ Z^A^ for our e+e~ Z^A^ process. Although hdecay 
does not have this final state for the pseudoscalar Higgs for comparison, we did calculate the branching ratio Ti{A^ — > 
Z^Z^) and added it to our analysis. The matrix elements squared, listed below, have the opposite effect from the 
^Z^ channel as they become more important as the mass of the pseudoscalar grows. The matrix elements grow like 
positive powers of (M^o — 2Mz)"' ■ This can be seen in the positive powers of this mass difference in the first and 
third terms of the matrix elements. 
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SUSY 1 2 



1 



327r2 (167r2)2 



2 + 



AO 



2A/|)2 



+ Af4 



4A/| 



E^l 

E^l 



^^^^^^ - ^^l^)2Keal(E AKi^ir)}. 



(7) 



The results of these partial widths can be seen in Fig. These graphs were created using the output of hdecay 
and adjusted to include the new contributions from the squark loops in the A^ — > 7Z° channel as well as the entirely 
new A" Z^Z^ channel. 
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FIG. 2: Improved branching ratios for the pseudoscalar Higgs boson (including the Z^^ Z''^ channel). The top and bottom quark 
loops as well as stop and sbottom squark loops have been included in the yl" 7^° and A° Z°Z° channels. The other 
branching ratios have been taken from hdecay and have been adjusted to allow for the new channel. We see that the 66 channel 
dominates at large tan/3, but the Z'^Z^ channel also plays an important role. The dip in the Z^ Z^ channel is not a kinematic 
one, nor is it due to the opening of a new decay channel in this case (the W^W' channel has not been included). It is due to 
a value of the pseudoscalar Higgs mass that allows for large cancellations to occur. 





FIG. 3: Diagrams contributing to the e'^e~ — > Z'^A" process. The graphs with the top and bottom quarks in the loop and 
the box graph are referred to as the standard model contributions in the text. The squark loop graphs are the dominant 
contribution to the processes. 
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Using our expressions for the two three-point functions with squark loops, the squark contributions to the matrix 
elements for the e+e" — > Z'^A'^ process can be written 

, 3„3 .,2 (t^ ~MlM%) r / 2sMl 

Ml U {tu-MlM%:, 



'{t + uf 



where, 
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-Af|Mloj(^|x2|2 + |x4|' 

2 Real ^xiXj + X2,xX ]\, (8) 



2^2 



2^3 



2:4 



V^ + ^i^A' (9) 
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The diagrams contributing to the e'^e~ — > Z'^A'^ process are shown in Fig. (O. The process was written as 
e+(pi)e^(p2) ^ ^°(^P3)^"(— Ps) and we have employed the usual kinematic variables s = (pi +P2Y , t = {P2 +P3)^, 
and u = {p3 + pi)^. We took the electrons to be masslcss and set pf = P2 = 0, p| = , and p^ = AI^q. 

The differential cross-sections for a future ILC are shown in Fig. The squark contributions can be seen for the 
two different center-of-mass energies (500 GeV and 1 TeV). The quark contributions are not shown because they are 
several orders of magnitude smaller than the squark contributions which are on the order of 0.2 fb for M^o = 120 GeV 
and tan/3 = 1. The contributions from the SM fields and the MSSM fields would not have mixed in any of the processes 
in this paper due to their tensor structure, as described in the appendix. The differential cross-section is symmetric 
with respect to the scattering angle cos 9 and with the interchange of the kinematic variables t ^ u. 



V. CONCLUSIONS 



In conclusion, we have found that the dominant contributions to the e+e" Z^A^ process does not come from 
top and bottom quark loops. The extreme enhancement of the pseudoscalar Higgs-squark-squark vertices due to the 
trilinear soft-supersymmetry breaking terms make the squark loops the dominant contribution to this process. This 
result is generically true in the entire neutral Higgs sector of the MSSM. Based on these results, squark loops appear 
to play a very large roll in Higgs phenomenology in the neutral sector. 
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APPENDIX A: MATRIX ELEMENTS 



The matrix elements for the process e+e^ — > Z'^A^ has implicit in it two one-loop three-point functions '^Z^A^ and 
Z^Z'^A^. We constructed these two three-point functions with quark and squark contributions for the collider process 
and to understand the squark contributions to our decay widths. There are also four box-type quark diagrams that 
have not been discussed so far and are not related to the three-point functions analyzed in this paper. These were 
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FIG. 4: Differential cross-sections for different parameters and center-of-mass energies at a future international linear collider 
for the process e^e~ Z'^A'^. All are for the case of maximal mixing in each of the squark sectors (sin(26i,i,) = 1) which leads 
to the largest contributions. Only the contributions from the squark loops are shown as they are several orders of magnitude 
greater than the contributions from the SM fields which are on the order of 0.2 fb for Af^o = 120 GeV and tan/3 = 1. Increasing 
the mass of the pseudoscalar drops the differential cross-section by about 30% and increasing the center-of-mass energy to 1 TeV 
decreases the differential cross-section by about 75% for tan/3 = 1. 



included in our analysis but they are not part of the squark sectors for which explicit results will be presented and 
are negligible in size to the SM contributions. 

The matrix elements for the rare pseudoscalar decays A'^ ^Z^ and Z^Z^ can be written as the sum of two 

parts. These are one loop decays that can have standard model fields and supcrsymmctric fields in the loops. The 
standard model and supcrsymmctric contributions do not mix due to their tensor structure. The supcrsymmctric 
diagrams can be broken into two form factors based on gauge invariancc, for the A^ijp^) — > ^{—Q^)Z^{—p'Q decay we 
can write the one-loop ^Z^A^ vertex with squark loops as 

The Z'^Z^A^ one-loop vertex will take the same form, but will have different coefficients ^'^id ^V)- There is 
no term that is proportional to the e tensor in the squark loops because the pseudoscalar docs not couple to squarks 
with a 7'"'. However, this is the case when there are standard model fields in the loop. Therefore, the same ^Z'^A^ 
vertex with quark loops can be written as 

*ri;:sM = e^''"''p3,aQ/3. (A2) 

When we try to interfere these two vertices it is now easy to see that we find zero due to a repeated index in the e 
tensor in the first term and a term anti-symmetric in the indices multiplied by one symmetric in the indices in the 
second term. Therefore, the decay widths will be additive for these processes. This is also true for the differential 
cross-section. Thus we can write 

TtOT = r.SM + TsUSYj (-A3) 

c^CTtot = dusM + da svsY- (A4) 
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This leads to some interesting phenomenology. It needs to be noted that the pseudoscalar Higgs boson does not 
appear in the Standard Model, thus the standard model contributions listed here are the contributions from the 
standard model fields, in this case, the top and bottom quarks. Thus deviations from the SM contributions do not 
tell us about the existence of supersymmetry in nature, but they do tell us about the mixing in the squark sector 
directly. 

When we add all the squark loop diagrams we can determine the unknown coefficients in our vertices. Thus, we 
can write 



= l!^AMi/,2j,3[^^^(112) + C24(221) - 2C24(111) 



- 2Q,[i3o(21) - So(22)]] - Qlsl[Bo{22) - Bo{n)]j - ^ c,; 1 ^ 2 j, (A5) 
M = A^B„ ^ 0; C24 ^ C12 + C23), (A6) 



7 7 
for the 7Z"A° vertex and 



A| = ^!£j^|f^|c4(r3)2[C24(121) + C24(211)-2C24(lll)] 

- clT^Qqsl [C24(122) + C24(121) + ^24(211) + ^24(221) - 4C24(111)] 

+ 4c^T3[Bo(12) - Bo(22)] - 2Q,sl[Bo{22) - Bo(ll)] 

+ QlsiC2A{lll) - slcl{T^f [C24(122) + ^24(121) + ^24(112) + C24(221)] 



Sq^ Cq]l^2y (A7) 

El - A%{Bo ^ 0; C24 ^ C12 + C23) (A8) 

for the Z'^Z'^A^ vertex. The Cij functions are the usual functions that appear in the Passarino-Veltman reduction 
prescription ■ It is easy to see that these expressions are finite. The B^ functions appear in pairs with opposite 
signs, canceling the 1/e poles which are independent of their arguments. The C24 functions also have argument 
independent 1/e poles that all cancel in the expressions. The same is true of the C12 and C13 functions. These 
functions can be written out fully as 



C,,(123) = C^j{MlQ^Mlo;ml , ml,ml) (A9) 
2 2 ^ 

^qi ' '^92 ' 



Bo(12)^i?o(M|;m? ,m? ) (AlO) 



where ~ {0, A/|, s} for the A" 7^", A" Z^Z^ , and e+e Z'^ processes respectively. This completes the 
missing squark contributions to the decay widths and the differential cross-section for the e+e^ — > Z'^ AP process. 
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